function [A,B] = gammaAB(mu,sig)
% function [A,B] = gammaAB(mu,sig)
% Calculates the A and B parameters of a gamma distribution with mean mu
% and standard deviation sig

A = mu/((sig^2)/mu);
B = (sig^2)/mu;